ANY ORDER SUPERCONVERGENCE FINITE VOLUME SCHEMES 
FOR ID GENERAL ELLIPTIC EQUATIONS 

WAIXIANG CAO *, ZHIMIN ZHANG t, AND QINGSONG ZOU * 

Abstract. We present and analyze a finite volume scheme of arbitrary order for elliptic equations 
in the one-dimensional setting. In this scheme, the control volumes are constructed by using the Gauss 
points in subintervals of the underlying mesh. We provide a unified proof for the inf-sup condition, 
and show that our finite volume scheme has optimal convergence rate under the energy and L 2 norms 
of the approximate error. Furthermore, we prove that the derivative error is superconvergent at all 
Gauss points and in some special case, the convergence rate can reach h 2r , where r is the polynomial 
degree of the trial space. All theoretical results are justified by numerical tests. 

1. Introduction. The finite volume method (FVM) attracted a lot of attentions 
during the past several decades, we refer to[l[3[i[71[TTJ[I71[ia[Tl[2Il[31[23[lHl[31] 

and the references cited therein for an incomplete list of references. Due to the local 
conservation of numerical fluxes, the capability to deal with the problems on the 
domains with complex geometries, and other advantages, FVM has a wide range of 
applications in scientific and engineering computations (see, e.g., [IS ! l2~T ] ). 

There have been many studies of the mathematical theory for FVM, see, e.g., 
[U [35] and the monographs [5J [TH1 US]. However, much attention has been paid to 
linear FVM schemes(see e.g., [S HI ED 1231 Hi]), high order FVM schemes have not 
been investigated as much or as satisfactorily as linear FVM schemes. Moreover, 
the analysis of high order FVM schemes are often done case by case. For instances, 
earlier works on quadratic FVM schemes can be traced back to [25j [16j ED], high 
order FVMs for ID elliptic equations were derived in [22 , and high order FVMs on 
rectangular meshes were derived and analyzed in [7], the quadratic FVM schemes on 
triangular meshes have also been intensively studied by [19j [28j [10] . To the best of 
our knowledge, no analysis of FVM scheme of an arbitrary order has been published 
yet. 

In this paper, we study a family of any order FVM schemes in the one-dimensional 
setting. Instead of a case-by-case study as in the literature for quadratic and cubic 
FVM schemes, we adopt a unified approach to establish the inf-sup condition. Earlier 
works(see, e.g. [20j [19] [28j [TTJ ) utilized element-wise analysis to prove that the bilinear 
form resulting from FVM is positive definite, which is a stronger property than the 
inf-sup condition. Hence, some assumption is needed for the mesh, such as quasi- 
uniformity and shape-regularity (in 2D). The major difference here is that we prove 
only the inf-sup condition (instead of positive definiteness of the bilinear form) and 
there is no mesh condition attached. With help of the inf-sup condition, we obtain 
the optimal rate of convergence in both the H 1 and L 2 norms. 

In this paper, we also study the superconvergence property of our FVM schemes. 
Note that while the superconvergence theory of finite element methods (FEM) has 
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reached its maturity ([21111126, 32 , 35]), the superconvergence analysis of FVM has also 
been focused on the linear schemes (see, e.g., [3 [55]). in this paper, it is shown that 
for a ID general elliptic equation, the superconvergence behavior of FVM is similar to 
that of the counterpart finite element method. For instances, the convergence rate at 
nodal points is h 2r , the rate at interior Lobatto points (defined in Section 4) is h r+2 , 
the convergence rate of the derivative error at Gauss points is h r+l . While in some 
special cases, some surprising superconvergence results have been found and proved. 
That is, the convergence rate of the derivative error at all Gauss points can reach 
lf+2 Qr ^2r ^ depending on the coefficient of the elliptic equations. The derivative 
convergence rate h 2r doubles the global optimal rate h r , which is much better than 
the counterpart finite element method's h r+1 rate; the derivative convergence rate 
h r+2 is one order higher than the counterpart finite element method's h r+1 . 

We organize the rest of the paper as follows. In Section 2 we present an arbitrary 
order FVM scheme for elliptic equations in one-dimensional setting. In particular, we 
use the Gauss points of order r > 1 to construct the control volumes and choose the 
trial space as the Lagrange finite element of rth order with the interpolation points 
being the Lobatto points of order r. In Section 3 we provide a unified proof for the 
inf-sup condition and establish the optimal convergence rate both under H 1 and L 2 
norms. In Section 4, we study the superconvergence property at some special points 
of our FVM schemes of any order. In Section 5, a post-processing technique based 
on superconvergence results in the section 4 is proposed to recover the derivative of 
the solution. Numerical experiments supporting our theory are presented in Section 
6. Some concluding remarks are provided in Section 7. 

In the rest of this paper, li A < B" means that A can be bounded by B multiplied 
by a constant which is independent of the parameters which A and B may depend 
on. 11 A ~ S" means "A < B" and "B < A" . 



2. FVM schemes of any order. In this section, we develop finite volume 
schemes for the following two-point boundary value problem on the interval f2 = (a, b): 

-(au')'(x) + p(x)u'(x) + y{x)u(x) = f(x), Vx G f2, 

u(a) = u(b) =0, ' ' 

where a > ao > 0, 7 — |/3' > k > 0, a, f3, 7 G L°°(fi), / is a real-valued function 
defined on Ct. 

We first introduce the primal partition and its corresponding trial space. For a 
positive integer N, let Zjy := {1, • • • , N} and a ~ xo < xi < . . . < xn = b be N + 1 
distinct points on Q. For all i £ Zjv, we denote t; = [xi-\,Xi] and hi — Xi — Xi-i, let 
h = max hi and 

V = {n : i e Zn} 

be a partition of 57. The corresponding trial space is chosen as the Lagrange finite 
element of rth order, r > 1, defined by 

U , v = {ve C(tt) : v\ T . G P r ,Vrj G V,v\ dn - 0}, 

where P r is the set of all polynomials of degree no more than r. Obviously, dim Ui — 
Nr — 1. 
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We next present a dual partition and its corresponding test space. Let G\, . . . , G r 
be r Gauss points, i.e., zeros of the Legendre polynomial of rth degree, on the interval 
[—1,1]. The Gauss points on each interval n are defined as the affine transformations 
of Gj to n, that is, 

9i,j = ^(xi + Xi-i + hGj), j € Z r . 
With these Gauss points, we construct a dual partition 

V = K ,r^ r } U {-rij : (i,j) eZ N x Z r J, 



where 



here 



T i,0 = [°)fl , l,l]' T W,r = [ffJV.r, l],Tj = [#i,j, fft,j+l] 5 



r if i e Z;v-l j w . ™ 

r-1 if i = TV ffi.r+i =ffi+i,i,V« S Aiv-i. 

The test space Vpi consists of the piecewise constant functions with respect to the 
partition V 1 , which vanish on the intervals r{ U t' n . In other words, 

Vpi = SpanjV'i.j : (i, j) S Zjv X Z ri } , 

where = X[g t j- + i] is the characteristic function on the interval r/j. We find 
that dimVp/ = Nr-1 = dim U£. 

We are ready to present our finite volume schemes. Integrating (|2.1[) on each 
control volume [gij, 9i,j+i], (hi) S Zjv x Z ri yields 

—(au')'(x) + (3(x)u'(x) + j(x)u(x)dx = / f(x)dx. 

9i.j ->9i,i 

In other words, 

«(9ij)"'(3ij)- Q! (3ij+i) u '(5ij+i)+ / (/3(^)it'(x)+7(x)u(a;))drc = / f(x)dx. 

(2.2) 

Let wp< S Vp' , wp/ can be represented as 

iV r< 
i=l j=l 

where w[ jS are constants. Multiplying (|2.2[) with tOjj and then summing up for all 
G Zjy x Z r( , we obtain 

2_, Z_, "''•< a (9i,i) u '(9i,j) ~ a(9i,j+iW(gi,j+i) + / ((3(x)u'(x) + ~f(x)u(x))dx 

i=l j=l \ J 9i,j J 

b 

f(x)w-p> (x)dx, 



or equivalently, 

N r N r 



22z2[ Wi <A a (9ij) u '(9i,j) +22z2 Wl <i ( / (^(x)u'(x)+j(x)u(x))dx 

i=l j=l i=l 3 = 1 \ 9i '3 y 



f(x)w-p> (x)dx, 

where [1^,3] = — w i,j— l is the jump of w at the point gjj, € Z^v x Z r with 
wi.o = 0, Wn,t — and Wija = Wi—i,v> 2 < i < N . 

We define the FVM bilinear form for all v G (f2), w-p' G Vp' by 

JV r 



apt". «>7") = ' ir ' ;// '-^ 

»=1 3=1 

+ 2^2^^,j ly (^(«)i;'(ar)+7(a:)t;(a:))di 



i=l j=l \ J 9i,j 



The finite volume method for solving equation (|2.1|) reads as : Find u-p G U£, 
such that 

a-p(u-p,w-p>) = (/, w-pi), Vw-pi e Vp/. (2-4) 



3. Convergence Analysis. In this section, we prove the inf-sup condition and 
use it to establish some convergence properties of the FVM solution. 

3.1. Inf-sup condition. We begin with some notations and definitions. First 
we introduce the broken Sobolev space 

W™' P (Q) = {v G C(Sl) :v\ Ti G W m ' p (Q),y T i G V}, 

where m is a given positive integer and 1 < p < oo. When p = 2, we denote H^} for 
simplicity. For all j > 0, we define a semi- norm by 



E ML, 



and a norm by 



\v\\ m , P ,v = ( E l^'L.P 

We often use | • \ m j> instead of | • \ m ,i,v and |j • \\ m ,-p instead of || • ||m,2,p for simplicity. 

n n 

Secondly, for all w-pi G Vp>, w-pi = X) 2 w i,j' l l } i,3i ^ 

i=l3=l 

N r N n 

\ W -P'\l,V' = EE A « :1 [ Wi jl 2 ' ll W Hlo,7" = EE^ W L 

4 — 1 j—1 l — l j — 1 



and 

II w HIp' = l w Hi,-p< + IkHlo.-p'- 

Noticing that wi.o = wjv,r — 0, it is easy to obtain the following Poincare type 
inequality 

Ihp'llo.P' ~ l w Hi,-p" Vw ?" e v T»-> ( 3A ) 

where the hidden constant depends only on SI and r. 

Thirdly, we denote by Aj , j £ Z r the weights of the Gauss quadrature 

r 



for computing the integral 



/(F) = J F(x)da 



It is well-known that Q r (F) = 1(F) for all F £ P 2r _i(-1, 1). Naturally, the weights 
on interval Tj, i £ Zjv are 



Aij — — Aj, j £ Z r . 



Before the presentation of the inf-sup condition, we define a linear mapping lip : 
UJp^V v , by 

u v v v = ^2^2v it ji/ji t j, 

i=l j=l 

where the coefficients Vij are determined by the constraints 

[vij] = A itj v' P (gi ! j), £ Z N x Z ri . 

Note that v-p £ U^,, the derivative v' v £ P r _i(ri), i eZn, then 
N r „ b 

^2^2 Ai,jv'-p{gi,j) = / v' v (x)dx = v-p(b) - v-p(a) = 0. 

i=l j=l Ja 

Therefore, 

N n N r 

v N ,r-i = ^2^2[vi,j] =^2^2 A i,jv' v (gi,j) - A N . r v' r {g N . r ) 

i=l 3=1 i=l j=l 

= —AN, r Vp(gN,r)- 

In other words, we also have 

[vN,r] = VN,r ~ UjV.r- 1 = A N ^ r v' v (g N ^ r ). 
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Consequently, 

N r N r „fc 

inpMi,-P'=EEv 1 [ t, «] a = EEv 1 (^j^(ff«)) a ~ / {v'vWfdx. 

i=l 3=1 i=l j=l , ' a 

Namely, we have 

\U-pvv\i,v ~ |t>7»|i,7»< (3.2) 

With all these preparations, we are now ready to present the inf-sup condition of 
a v (-, •)■ 

Theorem 3.1. Assume that the mesh size h is sufficiently small, then 
. a v (v v ,w v ,) 

mi sup - — — > c , (3.3) 

»PE^» pl 6V \\vv\\l,v\\Wv'\\v 

where cq > is a constant depending only on r,ao, k and tt. 
Proof. It follows from the bilinear form (|2.3[) that 



a-p(v-p, U-pv-p) = h + L 2l Vup S Up 

with 

N r v^v^ [ 9i ' i+1 

h = ^^[tijjay^tffij), h = E E u m / ( ( 5(a;)up(a;)+7(a;)wp(x))da;. 

i=l j=l i=l j = l "^.J 

Since (v' P ) 2 £ P 2r -2(T,),ie Zat, we have 



p {v' v {x)fdx = j2A^{v' v {g^)Y 



Therefore, 

iV r 
i=l j=l 

We next estimate 1-2. Let F(x) = (fi{s)v'-p{s) + 7(s)i>-p(s)) (is and denote by 
Ei= I v' v (x)V(x)dx - E A ii jv' r (gij)V(g ii j), 
the error of Gauss quadrature in the interval Tj , i S Zjv ■ Then 

W r b N 

h = - EE V ^i,i) = - Mx)V(x)dx + J2 E i- 

i=l j=l i=l 

Using the fact that v-p(a) — v-p(b) = and 

b l nb 

f}(x)v' v {x)v-p{x)dx — —— / (3' (x)vp(x)dx, 

^ .In 



we obtain 



v'-p(x)V(x)dx = / I j(x) 



[3'{x) 



v^(x)dx > k\\v v \\1. (3.4) 



On the other hand, by[14](p98, (2.7.12)), for all i G Z N 

^ +1 (H) 4 



Ei = 



(2r + l)[(2r)!] 3 

where ^ S t^. By the Leibnitz formula of derivatives, we have 

2r 



fc=r+l ^ ' 



with 



ci =max{||^||| T ._ 1>0O)Tj ,||7||i r ._ 1)OOiT4 } ^ ( fc J- 



Therefore, by the inverse inequality that 

||fp||r,oo 



we have 



ci(r!) 



(2r + l)[(2r)!] 3 
Combining the above estimates, we obtain 

ci(H) 4 



h > k\\vv\\ ,v 



(2r + l)[(2r)!] 3 
Then for sufficiently small h, we have 



h 2 \v-p\i-p ■ 



/ tt \ CKq . .9 ^ii ii2 • r tii 1 1 2 

Noticing (|3.1[) and (|3.2j) . we obtain 

\\Ii v v v \\ v , < \\vr\\i,v 

Therefore, for any v-p S 

a-piv-p ,w-p>) a-p(v-p,H-pv-p) 



sup 

w v ,ev v , \\w V '\\p> 



\U v vp\ 



> co\\vv\\i,v, 



where cq is a constant depending only on r, cto,n and J7. The inf-sup condition (13.31) 
follows. □ 

Remark 3.2. In the above proof, the partition V does not need to satisfy any 
quasi-uniform or shape-regularity condition. Moreover, (|3.3[) holds even the order of 
polynomials in each sub-interval Ti are different. 
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3.2. Energy norm error estimate. Following j^H], we use the inf-sup condi- 
tion (|3.3p and the framework of Petrov-Galerkin method to present and analyze the 
finite volume element method (|2.4I) . 

We first introduce a semi- norm and a norm in the broken space iJp(Q) = W-p (ft) 

by 



Ml=J2 \<T,+h*\v\l^ \M 2 v = \\v\\lv + \v\v 



It is straightforward to show that, 

\vp\-p ~ \vv\i,v, \\vv\\v ~ \\vv\\i,T, Vf-p G U?p. 
With these equivalences, the inf-sup condition (|3.3j) can be written as 



. a v (v v ,w v >) 

mi sup j — — n — > c 2 , (3.5) 

vv£U!p Wv , eVv , \\vp\\-p\\wp>\\v' 

where c 2 > also depends only on r,ao,n and fl. Moreover, we define a discrete 
semi- norm | • | 1 for all v G i?o (^) by 



lG,l 



AT r 



We next discuss the relationship between | • |-p and | • \c,i- First 

|v-p|g,i = |"p|i,-p ~ \vv\v, V v v G C^-p- 
On the other hand, for all v G flo(fl) n ilp(fi), 

(«'G^)) 2 < V'llw'lli^n) + M»"[|£ a( rO. (m) e Z ^ x Z - 
where the hidden constant depends only on r. Thus by the fact Aij < hi, we have 

N r 

i=i j=i 

A? 



."II 2 



Namely, 



where the hidden constant depends only on r. 
We are ready to show our main result. 

Theorem 3.3. Assume that u is the solution of (|2.1|) . u-p is the solution of 
Then the finite volume bilinear form a-p(-, •) is variationally exact: 

a v (u,w-pt) = (f,wv>), V w-pi G Vp>, (3.6) 



and bounded : for all v £ Hq(Q) D iJp(fi), wp> G Vp> , 

\a v {v,w v >)\ < M\\v\\ v \\w v >\\ v >, (3.7) 
where the constant M > depends only on r, cvq, k and f2. Consequently, 

M 

u — u-p L, < — inf u — up (3.8) 

117 c 2 «7=e^p ' 



where c 2 is i/ie same as m 

Proof. First, the formula (13.61) follows by multiplying (|2.1I) with an arbitrary 
function u>-p< £ Vp> and then using Newton-Leibniz formula on each control volume 
[9i,j,9i,j+i], (hj) 6 2jv x Z n . 

Next we show (|3.7j) . By the Cauchy-Schwartz inequality, from (|2.3|) there holds 
for all v e H&(Sl) H H v (n), w v > e Vp< that 



J=l ,J 

l 

(iV r< \ 5 / JV \ 2 

EE''-"!. I (Di^ + ii^jJ 

where the constant M depends only on r, a , k and Q. 

Finally, combining the inf-sup condition (|3.5p . (|3.6p and (|3.7p . we derive 
following similar arguments as in Babuska and Aziz ([2]), or Xu and Zikatanov (|27j). 
□ 

Corollary 3.4. Assume that u e i?o(0) n IPp +1 (fl) is the solution of (12. ip . 
and u-p is the solution of FY M scheme (|2.4|) . then 



u v \\i <h r \u\ r+ i^ v , (3.9) 

where the hidden constant is independent of "P. 

Proof. It follows from the definition of || ■ \\p and (|3.8p that 

\\u — Uj>\\i < \\u — up\\p < inf ||u — v-p\\-p. 

Noticing that 

/ N \i 
inf ||u - vp\\p < \\u - ui\\i + V] ft?|u - U/I2 r , 

where ui is the Lagrange interpolation of u at the Lobatto points (defined in the next 
section) in the trial space U-p. By the standard approximation theory, we obtain the 
estimate ([3T9"]> . □ 

9 



4. Superconvergence. In this section, we will present the superconvergence 
properties of the FVM solution. To this end, we need to use the interpolation of 
a function on the so-called Lobatto points. This kind of interpolation has been 
used in the literature for superconvergence analysis, see, e.g., [30l [31]. We denote 
Li, L2, • • • , L r -\ the zeros of Pr(x), where P r (x) is the Legendre polynomial of order 
r. Moreover, we denote Lq = — 1, L r = 1 and N r = {0, 1, • • • , r} for r > 1. The family 
of points Lj, j G N r are called Lobatto points of degree r. The Lobatto points on r, 
are defined as the affine transformations of Lj to Ti, i.e, 

k.j = ^(xi + Xi-i + hLj), j G N r . 
Let ui G Up be the interpolation of u such that 

ui(hj) = u(lij), eZjv x N r , 

then by [32](pl46, (1.2)) 

\(u - ui)'(gi,j)\ < //>lr+2,i,<<i' (4.1) 

where uj' 1} = ://,., i- //,.,<• 

Theorem 4.1. Let u G ^(fl) fl H 7 V +2 (Q) be the solution of (|2.ip . and uv the 
solution of FVM scheme (|2.4p . Then 

\a V (u - Ui,W-p')\ < h r+1 (\u\ r+2 ,V + \u\ r +l,oo,v) ||wp'|[t") ^ w V G Vp> • (4.2) 
Consequently, 

\\ U I ~ U"p\\l < h r+1 {\u\ r+2 .v + |u| r +i,oo,p) • (4-3) 

Proof. Recall the definition of bilinear form (|2.3|) , integral by part we obtain 

a-p(u - ui, w-pi) = h + I2 

with 

N r 

i=i 3=1 

/ 2 = EE^/ (7 -/?')(« -«/)• 

i I 3=1 

For alH G Zjv, note that 

+1 1 
\v. — ui){g i ^)<h r i + \u\ r+lyOCJyTz , \u\ r+2AyTt < h?\u\ r+ 2, n , 

by the Cauchy- Schwartz inequality, (14.11) and the standard approximation theory, 

(N r 
EE (^ (r+1 Vlr + 2,2 )Ti + ^ r+3 kl? + l,oo,r 4 + ll« " «/|lo,r 4 ) 
i=l 3=1 

< h r+1 (\u\ r +2,V + |«|r+l,oo,7') II^'IIt"- 
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The desired result (|4.2p is proved. 

We next show (|4.3p . By the inf-sup condition (13. 3[) . 

a-p(u-p — ui,w-p>) a-ptu — ui,w-p>) 

CO «/-itp i< sup n (i = sup n n . 

w v ,zv v , Wp'Wv w v ,ev v , Wwv'Wv 

Combining the above inequality with (|4.2|) . we derive (|4.3p . □ 



As a direct consequence of (|4.3|) . we have the following L 2 error estimate. 
Corollary 4.2. Let u <E #o(^) nifp +2 (fi) be the solution of (|2.ip , and up the 
solution of FVM scheme (12.41) . i/ien 

h-^||0<^ +1 |k||r+2,P, (4.4) 

where the hidden positive constant is independent of V . 
Proof. By the triangle inequality, 

II" - wp||o < ||u - w/||o + - wp||o- 
By the Poincare inequality and (I4.3p . we have 

- w^Ho ^ |w/ - < h r+1 \\u\\ r+2 ,v- 
Moreover, by the standard approximation theory, 

ll« - «/||o < h r+1 \\u\\ r+i < h r+1 \\u\\ r+2 , v . 

The desired estimate (|4.4p follows. □ 

Remark 4.3. In the above L 2 error estimate, we do not use the so-called Aubin- 
Nitsche technique. However, we need slightly stronger regularity assumption on the 
exact solution u. 

We next study the superconvergence property at the nodes Xi,i G Zjv_i. 
Theorem 4.4. Let u be the solution of (|2.ip . and u-p the solution of FVM scheme 
dI21) . Ifu G W%" +1,OD {Sl), ften 

2r+l 

\{u - up)(xi)\ < h 2r ^2 \ u \],oo,v, Vi G Zjv_i. (4.5) 

j=r+l 

Proof. Let e = u — u-p and 

e(a:) = / {P(y)e'(y) + j(y)e(y))dy, Vx e [a, 6]. 



By the construction of the FVM scheme, both u and u-p satisfy (12.21) . then for all 

G x Z n , 

-(«(5«+i)e'fe,j+i) - a(9i,j)e'(5<,i)) + - = 0. 

Namely, 

a{9i,j)e'{9i,j) ~ = Co, (4.6) 

where Co is a constant independent of i,j. 
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On the other hand, let G(-, •) be the Green function for the problem (|2.1[) . Then 
for all v G H^(n), 

a{y)v\y) — (x,y)dy + J (/3(y)v' \y) + j(y)v(y))G(x,y)dy. 
In particular, for all i G Zjv-i, 



e(xi) 



dG 

a(y)e'(y) — (xi,y)dy+ \ (/3(y)e'(y) + j(y)e(y))G(xi, y)dy. 



Noting that G{x l ,a) = G(x u b) = 0, byflU]) 
e{xi) -- 



dG 



a 

N r 



(y)e'(y) - e{y)) — {xi,y)dy 



dG 



E Ak 'i ( a (9k,j)e'(g k j) ~ s{jg k j))—{x u g kd ) + E x 
fc=ij=i v 

/b QQ 
— {x l7 y)dy + E 1 +E 2 =E 1 +E 2 , 



where 



1 ^ (2r+ l)[(2r)!] 3 



fc=i 



dG 

((a(y)e'(y)-s(y)) — (xi,y) 



(2r) 



^(2r + l)[(2r)!]3 



9G 
<9y 



l(2r) 



with ^ fc ,r/ fc G r fe . 

We next estimate E\ and £2, respectively. Note that &$> — «w for j > r and 
the Green function G(xi, •) has bounded derivatives of any order on each r k , k G Zjy, 
then 



AT 



2r+l 



fc=l V j— j— r-f 

iV / r 

;$E^ r+1 E fc_i ii e iio.~.n + E 

k=l \j=0 j=: 

2r+l 

<h T \\e\\ ^ v + h 2r+1 E K\oo,7> 



2r+l 



j=r+l 



j=r+l 

where in the second inequality we have used the fact that[l 

\e\j,oo,r k < /ifc J ||e||o,oo,r fc + /ife +1 "- ? |e|r+i,oo,r ii! Vj G Z r . 

We next consider the term ||e||o.oo.p- By the triangular inequality and the inverse 
inequality 

\\ui - ur\\o,oo,v ^ h?\ui - up\i < h*(\u - + \ui - U V \i) 

< h r+ ^\u\ r+1 < /i r + ^|ti| r +l,oo,P, 
12 



wc have 



||e||o,oo,-p < II" - u/||o,oo,7> + II"/ - Up\\o,oo,p < h r+ 2 \u\ r+ i 



oo.V- 



Therefore, 



2r+l 
j=r+l 

As for E2, a direct calculation shows that 

\E 2 \ <h 2r \\G\\ 2rt00 , v . 

Combining E\ with E2, we obtain (|4.5I) . □ 
As a direct consequence of (|4.5|) . we have 



E 



node 



(4.7) 



We next estimate the term e(xt) — e(xi-i),i G In which plays a critical role in 
our later superconvergence analysis. 
Theorem 4.5. For all i e 7L N , 



2r+l 



|(ti - u v ){xi) - (u - uv)(xi-i)\ < h 2r+1 22 

j=r + l 



(4.8) 



Proof. By the same arguments as Theorem 14. 4[ we obtain 



N 



e(xi) - e(xi-i) = ^2(E' hk + E' 2>k ), 



k=l 



where 



E i,k - 



E 2,k = 



h 2 k r+1 (rl) 



((a(y)e'(y) - e(y)) (j^( x i>y) ~ ~^j^ Xi -^ v '\ 



(2r + l)[(2r)!] 3 

h 2 k r+1 (r\) 4 (dG ,^_ ( n 
(2r+l)[(2r)!] 3 \dy [Xi,V) dy (Xl ~ 1,y) 



(2r) 



(2r) 



with ^ fc ,77 fc e r fc . 

Recall the construction of the Green function G(xt, •), for all j G N27., 

<9G dG 
\\-g-(xi,y) - -Q^(xi-i,y)\\j,oo,Q\Ti ^ h\\G\\ j+ i t00 ,n\ Ti , 



and 



dG 



ay 
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Since the Green function G{xi, ■) G C (ifc), k G Z N is bounded, then 



((a(y)e'(y) - e(u)) f^(xi,y) - ^(x t - U y) 



(2r; 



y=Zk 



2r 
3=0 



,dG, 



\ae — s 



dG, 
dy 



< 



>* E 



2r+l 



E 



VJ=0 



j=r+l 



Following the same estimate for X^=o l e li,°o,Tj» \ u \],oo,r k 
obtain 

which yields the inequality (|4.8p directly. □ 



as Theorem l4.41 we 



Next we present the superconvergence property of u'-p at Gauss points, and u-p at 
Lobatto points. Before our analysis, we first introduce a special polynomial. For all 
v[t) G ff x ([-l, 1]), we denote by 

r 

Vr (t) = ^2b j M j (t) 
3=0 

the rth approximation of v(t) with 

v(l) +v(-l) _ v(l)-v(-l) 



-, h = 



1 



bj = U-n) / v'm^^dt, J =2, 



where Mj is the Lobatto polynomial of degree i and Lj is the Legendre polynomial of 
degree j . For all x G , i G Z jv , we denote by 

Vr (*) = « r (^±^I±^), te[-l,l] 

the rth approximation of v(x) on the interval Ti,i G Zjy. Then 

\(V - V r )(h,j)\ < h r+2 \\v\\ r+2 ^,V, ieZr-l, (4.9) 

and 

|(«-«r)'(5ij)l S^NIr+^oo,*, iSZ r . (4.10) 

Theorem 4.6. Lei u G W£ +2: °°(f7) &e tfie solution of (|2~Tj) , and tip i/ie solution 
of FVM scheme flU}. TTien 

(4.11) 



|(u - M-p)'(5ij)| < /i r+1 ||u|| r +2,oo,P, G Z W X Z r _i 
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and 



\(u - u r ){l M )\ < h r+2 \\u\\ r+2 ^ v , e Z w x Z r . (4.12) 

Proof. For all v G H&(Sl), let 

A(u, v) — (cm', ?/) + (/W + ju, v). 

Then we have 

v(x) = A(v,G(x,-)), Viefi, 

where G(x, •) is the Green function for the problem (|2.1j) . Let Gp be the Garlerkin 
approximation of G(x, •), that is 

vp{x) — A{vp,Gp), \fv-p € C/p,Vx € Q. 

Thcn(see[9](p33)) 

- U r , Gp) < /i r+2 ||u||r+2,cx),7'- 

We next estimate the term A(u — up, G-p). Note that Gp £ UJp, then 
A{u-up,G v )= f (a(y)e'{y)~e(y))G' v {y)dy 

J a 

N r 



^2^2 A k,j( a {9k,j)e'(gk,j) -e(g k ,j))G v (g k)j )+E 3 
k=ij=i 

b 

3, 



= C [ G' v dy + E 3 = E : 

J a 

where e(y),s(y) and Co are the same as in Theorem 14. 41 and 

^[(a(y)e'(y)-e(y))G' v (y)] 



N i2r+l/ i\4 



^(2r+l)[(2r)!P 



y=ik 



Then 



JY 



£-3^^^fe r+1 I |Cylr,oo,-nJ e llr+2,oo,T fc + /] \@V \j,oo,T h l|e||g» — j*+2,< 
fc=l \ 3=1 

N / r-1 

<^2\Gp\ 2 ,i, Th hl +2 \\e\\ r+ 2,c a ,T k +^hl r+2 ~ J \\e\\ 2 r-j+2,oc,T k 



k=l \ j=l 

< /i r+2 ||e|| r+2 ,oo,P < h r+2 \\u\\ r+2 ^p. 
Here we have used (??) and the inverse inequality 

\Gp\j, 00 , Tk < /lfe _J |Gp| 2 ,l,r fc , Vj e Z r 
15 



and the fact [5] (p33) 



N 



\Gv\2,i,v — ^ \Gv\2,l,T k < C 



i=l 

with C a bounded Constance. Note that 

(u r — up){x) — A(u r — u-p, G-p) = A(u — u r , G-p) + A(u — up, G-p), 
then we have 

(u r - up)(x) < /i r+2 |M| r+ 2,oo,p. 

By inverse inequality, 

(u r - Up)'(x) < /i r+1 ||M|| r +2,oo,P- 

Combing above estimates with (|4.9p and (|4.10p . we obtain (14. 1 1[) and (|4.12[) directly 

by the triangular inequality. 

□ 

As a direct consequence , we have 



and 



where 



\u-up\a., <h r+1 , \u-up\ aver , 1 <h r+1 (4.13) 
\u-u v \ Lfl <h r+ \ \u-up\ aver , <h r+2 , (4.14) 



/ 1 N r 
Haver,! = ( J^YY^^j? 



i=l j=l 



and 



. N r \ z I N r 

\ V \ L fi = ( EE^^'n ' l U Uver,0= [jfr;YY V ( li >tf 
v »=l 0=0 J \ i=l j=0 

here Wij are weights of the Lobatto quadrature. 

Now we consider a special case that j3 = 0, we have the following theorem. 
Theorem 4.7. Let u be the solution of (|2.1jl . andup the solution of FVM scheme 
(H3D- IfP{x) = 0,Va; W^ +2 '°°(ft), t/i. 



en 

2r+l 

\(u-up)'( 9lJ )\<h mi < r+2 ^ \ u koo,v- (4.15) 

fc=r+l 

Proof. First, both u and up satisfy (|2.2[) . there holds for all (i,j) € x Z r _i 

that 

/•S«,3+i 

a(3i,i+i)e'(fl , i,j'+i) - a(9i,j)e'(9i,j) = / 7(x)e(a;)dx, 
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which yields 

e'(^ + i) = + -r^—, / 9 " J+1 j(x)e(x)dx. (4.16) 

On the other hand, 

e{xi) - e(xi-i) = / e'(y)dy = ^ A itj e'(gi tj ) + E l , 



whereby QI](p98, (2.7.12)), 



^ = (27TT)pW (e ' )(2, ' )(6) ~ ^l"'*^-™ 6 e ^ 

By Theorem (|4~5l and (|4.16|l . we obtain 

hie'(g it i) + hi / 7(x)e(x)da; < /\ 2r+1 ^ Mfc.oo.p- 

•> Xi -1 i i -i 



Then 

2r+l 



e '(ffij) / |7(x)e(x)|dx + /i 2r |w|fe,oo,p, Vj e Z r . 



We next estimate the term J^* ^ |7(x)e(a;)|da;. Note that 

Si, j / /•■/ , /-.'-• , \ 

7(x)e(x)| < UtIIoc I / \(u— ui)(x)\ dx + / \(uj - u v )(x)\ dx J 

where in the above inequalities we have used |u— u/| < /i r+1 HitHr+i.oo < /i r+1 |u|r-+2,oo,'P 
and \ui — u-p\ < \ui — u-p\i,-p < h r+ ? \u\ r+ 2^v- 
Therefore, 

2r+l 

\(u~u P y(g hJ )\<h™^ r + 2 ^ J2 Mw- 

k=r+l 

The proof is completed. □ 



In particular, when /3 = 7 = 0, we have a better result. 

Theorem 4.8. Let u e W^ +1,00 {Vl) be the solution of (|2T|) . and up G the 
solution of FVM scheme $2A$). If j3(x) = j(x) = Q,Vx £ ft, /or aZZ E Z N x Z r; 
we /ioue 

2r+l 

\ u '(9i,j) - u 'v(9i,j)\ S h 2r \u\k,oc,v- (4-17) 

fc=r + l 
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Proof. By (|4.6p . we denote the constant 

C = <x(9i,j)(u'v(9i,j) ~ U '(.9 4J ))- 
The fact M-p £ P r yields that 



e(a;i) - e(xi_i) = / e(i)cft 

T i-x r 

= ^2 A *,k e '(9i,k) + / u'(t)dt - y^A itk u'(g iik ). 

k = l Jx i-1 k=l 



By [T4](p98, (2.7.12)) 



/ u'(t)(ft-5^Ai,fcU^(ff* 

7 1 



' X *-l k=l 

and Theorem 14.51 we have 



< h 2r+1 \u\ 2r+1 . 



r 2r+l 
■j = l fc=r+l 

which yields (|4.17l) directly. □ 

Remark 4.9. We see that at the Gauss points, when j3 — 7 = 0, the derivative 
convergence rate h 2r doubles the global optimal rate h r , which is much better than the 
counterpart finite element method's h r+1 rate, when (3 = 0, the derivative convergence 
rate h r+2 is one order higher than the counterpart finite element method's h r+1 ; and 
at the nodal points, the convergence rate h 2r almost doubles the global optimal rate 
h r+1 and equals to the counterpart finite element method's h 2r rate; and at the Lobatto 
points, the convergence rate h r+2 is one order higher than the optimal global rate h r+1 , 
which is the same as the counterpart finite element method. 

5. Post processing. We observe from (|4.11|) . (|4.15j) and (|4.17|) that u' v approx- 
imates the derivative of the exact solution u pretty well at the Gauss points. In this 
subsection, we will recover u' in the whole domain SI. 

For alH = 1, . . . ,N — 1, we construct a function Vi G P 2r -i([xi-i, by letting 

Vt(gi,k) = u' v (gi tk ), l = i,i + l; k = 1,2, . . . ,r. 
Then we define for all x G Tj = [xj_i, Xi], i = 1, . . . , N, 

!vi(x), i = 1, 

\{vi{x) + Vi-i{x)), 2<i<N-l, 
v N -i(x), i = N. 

To study the approximation property of u, we note that in each [arj_i,£i+i], 

u'{x) = (L 2r -iu')(x) H ^-^ — Y[(x - g itj )(x - g i+ i tj ),£ G [xi-i,x i+ i] 

where the Lagrange interpolant 

i+l r 



x - 9v, r 



(L 2r -iu')(x) =^2^2u'(gij)wij(x),wi tj (x) = J| 
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Noting that 



l=i 3 = 1 

we have 

u'(x)-Vi(x) =X/2J( U ' - u 'v)(9i,j)wi,j(x) + — 7^T| — \\{x-9i,j){x-gi+\,j)- (5.1) 

l=i 3=1 ^ j = l 

Since for all i = i, i + 1, j = 1, . . . , r, we have 

\wi,j{x)\ < C r ,Vx £ [Xi-i,X i+ i], 

where c r is a constant depends only on r, we obtain by (I4.11j) . (|4. 15)) and (|4.17[) that 

2r+l 

\u'(x) - Vi(x)\ < h m 22 \u\k,oo,-p, 

where m = r + 1 for general elliptic equations, m = r + 2 if /3 = 0, and m — 2r if 
/? = 7 = 0. Consequently, we have 

2r+l 

\u'(x)-v(x)\ <<h m J2 k| fe ,oo,P,Vx£n. 

fc=r+l 

6. Numerical experiments. In this section, we present numerical examples to 
demonstrate the method and to verify the theoretical results proved in this paper. 

In our experiments, we solve the two-point boundary value problem (|2.1|) by the 
FVM scheme (|2.4|) with r — 4 or r = 5. The underlying meshes are obtained by 
subdividing 57 = (0, 1) to N = 2, 4, 8, 16, 32, 64 subintervals with equal sizes. 

Example 1. We consider the two-point boundary value problem (|2.ip with 

a(x) — e x , (3(x) = cosx, j(x) = x, Va; 6 fi, 
and / is chosen so that the exact solution of this problem is 

u(x) = sinx(a; 12 — a; 11 ). 

We list approximate errors under various (semi-) norms in Table 16.11 ( for the 
scheme r = 4 ) and Table 16.21 ( for the scheme r = 5 ) . 

To explicitly show the convergence rate of different approximate errors, we plot the 
error curves in Figures 16.11 and 16.21 We observe from Figure 16.11 that the convergence 
rate \u — u-p\\ is r and the convergence rate of \u — «-p||o is r + 1. In other words, the 
FVM approximate solution converges to the exact solution with optimal convergence 
rates under both for H 1 and L 2 norms, as predicted in (|3.9[) and (|4.4p . We also 
observe that the error \m — u-p\\ is of order r + 1, which confirms the convergence 
result in (|4.3p . The errors \u — u-p\ ave r,0i \u — u-p\l,o and E no d e are presented in Figure 
16.21 It is observed that \u — u-p\aver,o and \u — Up|z,,o converge with a degree r + 2 
which confirm the superconvergence property at Labatto points given in Theorem 14. 61 
Since E no d e converges with a rate 2r, it confirms our theory in Theorem 14.41 

Example 2. In this example, we test the convergence behavior of derivative error 
at Gauss points. We consider three cases of Equation (|2.ip . they are 
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Table 6.1 
r = 4 



N 


\\u - u v \\ 


\\u - U-p\\l 


\uj - U V \i 


\u - u v \ Lfi 


2 


1.8618e-03 


5.1201c-02 


5.2554e-03 


3.3420e-04 


4 


1.4386e-04 


7.2801e-03 


3.1271e-04 


9.8931e-06 


8 


5.9282e-06 


5.9099e-04 


1.1758e-05 


1.8624e-07 


16 


1.9882e-07 


3.9516e-05 


3.8485e-07 


3.0490e-09 


32 


6.3240e-09 


2.5119e-06 


1.2166e-08 


4.8197e-ll 


64 


1.9850e-10 


1.5766e-07 


3.8129e-10 


7.5536c-13 


N 


\u — U-p \aver,0 


\u - U-p\ G ,l 


\u — U-p\aver,l 




2 


2.1895e-04 


8.0770e-04 


5.4962e-04 


1.1874e-05 


4 


6.2680e-06 


5.3025c-05 


3.5877e-05 


5.9186e-08 


8 


1.1716c-07 


1.9692e-06 


1.3338e-06 


2.3666e-10 


16 


1.9150e-09 


6.3947e-08 


4.3328e-08 


9.2827e-13 


32 


3.0260c-ll 


2.0170e-09 


1.3667e-09 


— 


64 


4.7425e-13 


6.3175c-ll 


4.2809c-ll 


— 






lAbLb 


Z 








7* = 5 






N 


\\u - u v \\ 


\\u - U V \\l 


\UI - Up\l 


\u - u p \ Lfl 


2 


4.8206e-04 


1.5546c-02 


8.5017e-04 


4.0891e-05 


4 


1.5627e-05 


9.6503e-04 


2.1627e-05 


5.2643e-07 


8 


2.9713e-07 


3.6434e-05 


3.8065e-07 


4.6413e-09 


16 


4.8711e-09 


1.1927e-06 


6.1190e-09 


3.7318e-ll 


32 


7.7022e-ll 


3.7707e-08 


9.6282e-ll 


2.9365e-13 


64 


1.2073e-12 


1.1817e-09 


1.5081e-12 




N 


\u — Up \aver,0 


\u ~ u v \g,i 


\U — U-p\aver,l 


Enode 


2 


2.6075e-05 


2.2179e-04 


1.4819e-04 


4.6819e-08 


4 


3.2965c-07 


5.8493e-06 


3.9162e-06 


3.0508c-ll 


8 


2.8971e-09 


1.0085e-07 


6.7553e-08 


2.6318e-14 


16 


2.3277e-ll 


1.6089e-09 


1.0779&-09 




32 


1.8311e-13 


2.5266e-ll 


1.6928e-ll 




64 




3.9473e-13 


2.6482e-13 





Case 1 : a(x) = e x , j3{x) = cosx, j(x) = x; 
Case 2 : a(x) — e x , fi(x) = 0, 'f(x) = x; 
Case 3 : a{x) — e x , (3(x) — 0, 7(2:) = 0. 

The exact solution is always u(x) — sinx(a; 12 — a; 11 ) and the right-hand function / 
change according to the coefficients in different cases. 

Listed in Table 16.31 are errors in the derivative approximation at Gauss points 
for three different cases for r — 4 and r = 5, respectively. Plotted in Fig. 16.31 are 
corresponding error curves. We observe that the convergence rate is r + 1 for Case 1, 
r + 2 for Case 2 and 2r for Case 3. These numerical results are consistent with our 
theories derived in Section 4. 

7. Concluding remarks. The mathematical theory for the FVM has not been 
fully developed. The analysis in the literature for high-order FVM schemes are often 
done case by case. It is a challenging task to develop mathematical theory for FVM 
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Fig. 6.1. left: r = 4, right: r = 5 




Fig. 6.2. left: r = 4, rigfa: r = 5 

scheme of an arbitrary order. In this article, we provide a unified proof for the inf-sup 
condition of a family any order FVM schemes in one dimensional setting. Based on 
this, we show that the FVM solution converges to the exact solution with optimal 
order, both in H 1 and L 2 norm. 

We also studied the superconvergence of our FVM schemes. It is shown both 
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Table 6.3 
Gauss points. 





r = 4 


r — 5 


N 


Case 1 


Case 2 


Case 3 


Case 1 


Case 2 


Case 3 


1 


4.7633e-03 


4.1098e-03 


4.1493e-03 


1.5667e-03 


1.1105e-04 


1.0183e-04 


2 


5.4962e-04 


3.1457e-05 


2.9493e-05 


2.6075e-05 


1.9562e-06 


8.6701e-09 


4 


3.5877e-05 


4.2964e-07 


1.1316e-07 


3.9162e-06 


2.8751e-08 


3.1732e-ll 


8 


1.3338e-06 


8.4296e-09 


4.3677e-10 


6.7553e-08 


2.6812e-10 


3.6386e-14 


16 


4.3328e-08 


1.4052e-10 


1.7002c-12 


1.0779e-09 


2.1878e-12 




32 


1.3667e-09 


2.2319e-12 


6.6291e-15 


1.6928c-ll 


1.7284e-14 






Fig. 6.3. left: r = 4, right: r = 5 

theoretically and numerically that at the nodal and interior Lobatto points, the su- 
perconvergence behavior of FVM is similar to that of the counterpart finite element 
method. Moreover, in some special cases, the superconvergence property of the deriva- 
tive of the FVM solution at the Gauss points maybe much better than that of the 
counterpart finite element method. For instances, when /3 = 0, the convergence rate 
of the derivative of the FVM solution is h r+2 which is one order higher than the coun- 
terpart finite element method's h r+1 ; when f3 — 7 = 0, the order is h 2r which doubles 
the global optimal rate h r , and it is much better than the counterpart finite element 
method's h r+1 rate. In a recent study 29J, it is shown that after a simple post- 
processing procedure, the FEM solutions can have local conservation property. In 
this sense, the superconvergence property discovered in this paper become a powerful 
argument to support that the FVM still has its advantages. 
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